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ABSTRACT 

We examine how much information measured broad-hne widths add to virial BH mass estimates for 
flux hmited samples of quasars. We do this by comparing the BH masses estimates to those derived 
by randomly reassigning the quasar broad-line widths to different objects and re-calculating the BH 
mass. 

For 9000 BH masses derived from the H /3 line we find that the distributions of original and ran- 
domized BH masses in the MeH^redshift plane and the MBH~luminosity plane are formally identical. 
A 2D KS test does not find a difference at > 90% confidence. For the Mg II line (32000 quasars) we do 
find very significant differences between the randomized and original BH masses, but the amplitude 
of the difference is still small. The difference for the C IV line (14000 quasars) is 2 — 3cr and again the 
amplitude of the difference is small. Subdividing the data into redshift and luminosity bins we find 
that the median absolute difference in BH mass between the original and randomized data is 0.025, 
0.01 and 0.04 dex for H (3, Mg II and C IV respectively. The maximum absolute difference is always 
< 0.1 dex. 

We investigate whether our results are sensitive to corrections to Mg II virial masses, such as those 
suggested by Onken & Kollmeier (2008). These corrections do not influence our results, other than to 
reduce the significance of the difference between original and randomized BH masses to only 1 — 2cr 
for Mg II. Moreover, we demonstrate that the correlation between mass residuals and Eddington ratio 
discussed by Onken & Kollmeier are more directly attributable to the slope of the relation between 
H /3 and Mg II line width. 

The implication is that the measured quasar broad-line velocity widths provide little extra infor- 
mation, after allowing for the mean velocity width. In this case virial estimates are equivalent to 
Mbh oc L", with L/LEdd oc L^~°' (with a ~ 0.5). This leaves an unanswered question of why the 
accretion efficiency changes with luminosity in just the right way to keep the mean broad-line widths 
fixed as a function of luminosity. 

Subject headings: surveys — galaxies: active — quasars: emission lines — quasars: general 



1. INTRODUCTION 

It is now clear that accretion onto super-massive black 
holes (BHs) plays a crucial role in the process of galaxy 
formation and evolution. These BHs lie at the center of 
all massive galaxies and are thought to infiuence galaxy 
evolution via powerful radiative (i.e. quasar mode) and 
mechanical (radio mode) feedback (e.g. Croton et al. 
2006). The tight correlation between BH mass and the 
properties of the spheroidal component of their host 
galaxies (e.g. Tremaine et al. 2002) infers that the growth 
of both is closely related. 

The fundamental observables of black holes are mass 
and spin. While spin is still elusive (only being inferred 
from very indirect means), there are now various ap- 
proaches to the measurement of BH masses. The best 
constraint come from the BH at the centre of the Milky 
Way via stellar kinematics (e.g. Gillessen et al. 2009) and 
the BH in NGC 4258 via mega-maser orbits (Miyoshi et 
al. 1995). The masses of BHs in other local galaxies can 
be estimated by the impact of the BH on the motion of 
stars close to the nucleus, but model degeneracies remain 
(Gebhardt & Thomas 2009). 

For type-1 AGN, where we see all the way to the ac- 
tive nucleus, the high velocity gas surrounding the BH 
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gives US another probe of the BH potential. This is par- 
ticularly powerful because it does not rely on spatially 
resolving the sphere of influence of the BH and because 
the the light from a quasar often overwhelms the light 
emanating from its host galaxy. As a result BH masses 
can be estimated in high redshift quasars. The disadvan- 
tage to this approach is that gas kinematics can be more 
complex than the effectively collisionless stars. Physical 
processes such as radiation pressure and outflows (e.g. 
Marconi et al. 2008) may cause potential biases in BH 
mass determinations, although recent work suggests that 
this is not serious (Netzer & Marziani 2010). 

The approach to estimating BH masses using quasar 
broad emission lines relies on the assumption that the 
gas is largely in virial motion about the BH. If this is 
the case, we then require an estimate of the radius and 
velocity of the broad line region to estimate mass such 
that 
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where / is a geometric factor, -Rblr is the radius of the 
broad line region, (Tblp is the velocity width of the line 
and G is the gravitational constant. The only viable 
approach to measuring -Rblr is reverberation mapping 
(Blandford & McKee 1983). Evidence from multiple lines 
within the same objects is consistent with virial motion of 
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the broad line region (Peterson & Wandel 1999). Rever- 
beration mapping is limited to low redshift and low lumi- 
nosity broad line objects (i.e. Seyfert Is), largely because 
time-dilation and the scaling of -Rblr with luminosity 
(typically i?BLR ~ -^"'^ if driven by photo-ionization) 
can increase by an order of magnitude or more the time 
lag between continuum and line variability. In order to 
make estimates of BH mass in higher redshift quasars 
the correlation between -Rblr and luminosity found via 
reverberation mapping (e.g. Kaspi et al. 2000) has been 
applied. Using the i?BLR — L relation and calibrating to 
the H /3 emission line visible at low redshift, both the 
Mg II (McLure & Jarvis 2002) and C IV (Vestergaard 
et al. 2006) emission lines have been used for BH mass 
estimation. This approach, using the width of a broad 
optical or UV emission line and the i?BLR — L relation is 
generally termed the virial method. The power of such 
techniques is that they can be applied to large optical 
spectroscopic samples of quasars such as the Sloan Digi- 
tal Sky Survey (SDSS; Schneider et al. 2007), 2dF QSO 
Redshift Survey (2QZ; Croom et al. 2004) and 2dF-SDSS 
LRG and QSO Survey (2SLAQ; Croom et al. 2009). It 
is then in principle possible to study properties such as 
mass functions (e.g. Vestergaard et al. 2008) or the clus- 
tering of quasars as a function of BH mass (e.g. Shen et 
al. 2009). 

The above approaches give us great hope of charac- 
terizing the super-massive BH population over most of 
cosmic time. However, a number key issues have been 
raised. The C IV line in particular has been a subject 
of close scrutiny (e.g. Netzer et al. 2007; Baskin & Laor 
2005) largely due to the strong asymmetries and absorp- 
tion often seen in this line. However some authors argue 
that C IV can be used to provide a viable mass estimate 
(e.g. Vestergaard & Peterson 2006). This is because the 
asymmetries in C IV cause mass differences which are 
smaller than the relatively large scatter in the virial es- 
timates. Vestergaard & Peterson also argue that for the 
handful of objects which have reverberation mapping of 
multiple lines (including C IV), these provide consistent 
BH mass estimates. Not withstanding these points, the 
C IV line is thought to be located much closer to the 
nucleus than H /3 or Mg II so it is quite possible that 
different dynamics are in play in these two regimes. 

More generally, one of the recent challenges to the virial 
method is the study of the dispersion in quasar broad line 
widths as a function of luminosity for Mg II and C IV by 
Fine et al. (2008) and Fine et al. (2010) respectively. The 
most striking result from this work is the small scatter 
in line widths for the most luminous quasars, which is 
less than 0.1 dex for both the Mg II and C IV lines. This 
implies a scatter of less than 0.2 dex in BH mass. Given 
that this scatter must include dispersion due to accretion 
rate variations, orientation effects, scatter in the i?BLR — 
L relation and observational uncertainty, it is hard to see 
how such a small scatter can be produced. This scatter 
is also smaller than the typical quoted uncertainties on 
virial BH mass estimates. 

In this paper we consider the low scatter present in 
broad line velocity widths and address the question of 
how much information they actually provide concern- 
ing virial BH mass estimates. Our approach is to carry 
out a simple test of re-determining BH masses after ran- 
domly reassigning the broad-line velocity widths to dif- 
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Figure 1. The distribution of line width [log(D)], bolometric lu- 
minosity [log(LBoi)] and redshift for the samples analysed in this 
paper. The H /3 (black), Mg II (blue) and C IV (red) samples are 
shown separately. Contours are linearly spaced in 5 equal intervals. 
The lower left plot also shows the AGN with reverberation map- 
ping used to calibrate virial relations [open squares; taken from 
vestergaard & Peterson (2006)] and the approximate flux limits 
of the SDSS (top), 2QZ (middle) and 2SLAQ (bottom) using the 
mean bolometric correction from the B-band given by McLure & 
Dunlop (2004). For the 2QZ and 2SLAQ surveys the edge of the 
locus of sources lies close to these limits, while for SDSS which is 
i-band limited the sources have greater scatter across the nominal 
survey limits. Mg II and C IV line widths have been corrected 
to match the FWHM measurements of the H /? lines assuming the 
lines are single Gaussians. 

ferent quasars. Throughout this paper we assume a flat 
cosmology with Hq = 70kms^^Mpc~^, ilai — 0.3 and 
Qa = 0.7. 

2. BLACK HOLE MASSES 

The measurements of H /3 broad line widths is taken 
from Shen et al. (2008), who use the fifth data release 
of the SDSS quasar survey (Schneider et al 2007). Shen 
et al. largely follow the procedure of McLure & Dun- 
lop (2004), including their calibration of the virial BH 
mass estimator. This involves first fitting the underlying 
power-law continuum plus iron emission, which is then 
subtracted off. They then fit a double Gaussian to the 
emission line, constraining one component to be broad, 
and the other narrow. The FWHM of the line is then 
taken to be that of the broad Gaussian component. 

For the Mg II and C IV emission line fits we use the 
results of Fine et al. (2008) and Fine et al. (2010) re- 
spectively. The main reason for using the Fine et al. 
data sets is that as well as using quasars from the SDSS, 
they also use spectra from the 2dF QSO Redshift Sur- 
vey (2QZ; Croom et al. 2004) and the 2dF-SDSS LRG 
and QSO Survey (2SLAQ; Croom et al. 2009). This 
provides a larger dynamic range in luminosity at any 
given redshift. The line width measurements made by 
Fine et al. use the inter-percentile value (IPV) width 
(e.g. Whittle 1985). Continuum plus Fe II emission is 
first subtracted, and then a single Gaussian fitted. The 
IPV width is then measured over the range ±1.5 x the 
Gaussian FWHM. The IPV approach has the advantages 
that it is doesn't depend on a particular line shape, and 
as it is derived from the cumulative sum of flux in the 
line, easier quantification of the errors are possible (par- 
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Figure 2. a) The correlation between BH mass difference and Eddington ratio as given by OK08 using the data of Shen et al. (2008). b) 
The correlation between Mg II and H f) velocity width from the same data, which shows a gradient which deviates from 1-to-l (black solid 
line). The OLS(X|y) and OLS(y|X) fits are shown by the blue dotted lines and the bi-sector of these is given by the red dashed line, c) 
The relation between mass difference and Eddington ratio after correcting the Mg II velocity. 



ticularly incorporating the covariance due to continuum 
subtraction, which is generally not included in line fits). 
In the case of C IV Fine et al. (2010) also fit and sub- 
tract the contaminating He II and O III] at ~ I6OOA. 
The virial BH mass estimates use the McLure & Dunlop 
(2004) and Vestergaard & Peterson (2006) calibrations 
for Mg II and C IV respectively. The distribution of red- 
shift, bolometric luminosity and linewidth is shown is 
Fig. [TJ The fainter 2QZ and 2SLAQ samples allow us to 
probe a broader range of luminosities than just the SDSS 
alone, similar to the work of Kollmcicr et al. (2006) who 
also target faint AGN (hmitcd to i? < 21.5). 

The exact form of the virial relation and its calibration 
is not crucial in the analysis presented below. Our aim is 
to determine how much information the broad line veloc- 
ity width contributes to distribution of measured black 
hole masses. To this end we perform the simple test of 
randomizing the measured velocities. That is, assigning 
the velocity width of one quasar to another. We do this 
without reference to an object's redshift or luminosity 
(other than to only randomize H /3 widths with other 
H /3 widths, and similarly for Mg II and C IV), so that 
a quasar at the low luminosity and low redshift end of 
the sample could be given a width from a high redshift, 
high luminosity object. If the velocity widths are adding 
useful information to our black hole mass estimates, then 
the distribution of original black hole masses should be 
significantly different from the randomized black hole 



The work of Onken & KoUmeier (2008, OK08) sug- 
gested that Mg II based virial methods contain a bias 
which is a function of Eddington ratio. We also inves- 
tigate whether such a bias influences our results pre- 
sented below. OK08 present the bias as a mass differ- 
ence [M(Mg II)-M(H 13)] as a function of L/L^dd (Fig. 
[2^). However, these are both derived quantities which 
are correlated in non-trivial ways. In terms of measured 
(rather than derived) quantities this bias is due to the 
relation between H /3 and Mg II velocity width not hav- 
ing a gradient of 1 (see fig. [2b). OK08 define a matrix 
of corrections for the Mg II mass estimates. However, 
we find that a much more direct correction of the Mg II 
velocity widths also removes the bias. We find this cor- 
rection by making an ordinary least squares (OLS) fit to 
FWHM(Mg II) vs. FWHM (H /3) in both {X\Y) and 
(y|X) and taking the bisector (red dashed line in Fig. 
Wp; Isobe et al. 1990). The resulting fit is 

log[FWHM(MgH)] = 1.70 + 0.54 log[FWHM(H/3)]. (2) 

Applying this relation to correct the Mg II velocity 
widths results in an offset with respect to the original 
virial relation, which is corrected by using 
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= 0.505 + 0.62 log 
+2 log 
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1.278,(3) 



where the final term (+1.278) corrects the offset caused 
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Figure 3. The distribution of black hole mass estimates vs. redshift using the H f) (top left), Mg II (top right) and C IV (bottom left) 
lines. The black contours and points are the original BH masses (Shen et al. 2008 for H f} Fine et al., 2008 for Mg II and Fine et al. 2010 
for C IV), while the red contours and points with the broad line velocities randomized. The contours are used in regions of high density 
and are spaced at five equal linear intervals, using the same levels for both black and red contours. 



by the Mg II velocity correction. FWHM(MgII)c is the 
Mg II width corrected to match H (3 using Eq. [2l The 
result of applying the above BH mass estimate, including 
the correction to Mg II velocities, is shown in Fig. It 
can be seen that this corrects the strong bias originally 
presented by OK08. Applying such a correction as part 
of our analysis has no impact on the result presented in 
this paper. Further details will be described below. A 
gradient in FWHM(Mg II) vs. FWHM (H 13) which de- 
viates from one also has more broad implications for the 
virial technique. Assuming that the H /3 line is virialized 
and that the exponent for the velocity is 2 (as in Eq. 
[1]), then the implied exponent for the Mg II velocities is 
2/0.54 ~ 3.7. This deviation from a naive virial relation 
is similar to the results of Wang et al. (2009) who look at 
the calibration of virial relations without the assumption 
that the exponent for the velocity is 2. It could be that 
the virial assumption is wrong for either H /3 or Mg II or 
that the relative distribution of gas is markedly different. 
Infact, H /3 is expected to be due to multiple components 
(e.g. Sulentic et al. 2006). Our lack of understanding of 
the differences between H j3 and Mg II is a concern for 
the virial approach. 

3. RESULTS 

We first assess what the impact of randomization is 
on a single quasar by measuring the RMS of (Mq-Mr). 
We find RMSs of 0.54, 0.42 and 0.27 dex for H /3, Mg II 
and C IV respectively. In comparison to this the quoted 



Table 1 

The results of a 2-D KS test comparing original and randomized 
BH masses in both the Mbh ~ ^ a^nd Mbh ~ L planes. 



Line 


Test 








P{>D) 


H /3 (SDSS) 


Mbh - 


z 


8979 


0.0194 


2.18E-01 


H 13 (SDSS) 


Mbh - 


L 


8979 


0.0223 


9.38E-02 


Mg II (all) 


Mbh - 


z 


32214 


0.0290 


2.94E-08 


Mg II (all) 


Mbh - 


L 


32214 


0.0230 


8.27E-06 


Mg II (all, corrected)'' 


Mbh - 


z 


32214 


0.0280 


7.32E-02 


Mg II (aU, corrected)'' 


Mbh - 


L 


32214 


0.0236 


1.79E-01 


Mg II (SDSS) 


Mbh - 


z 


22910 


0.0295 


6.22E-06 


Mg II (SDSS) 


Mbh - 


L 


22910 


0.0259 


8.10E-05 


Mg II (2QZ) 


Mbh - 


z 


6784 


0.0221 


2.41E-01 


Mg II (2QZ) 


Mbh - 


L 


6784 


0.0212 


2.54E-01 


Mg II (2SLAQ) 


Mbh - 


z 


2491 


0.0252 


6.54E-01 


Mg II (2SLAQ) 


Mbh - 


L 


2491 


0.0278 


5.21E-01 


C IV (all) 


Mbh - 


z 


13795 


0.0213 


3.30E-02 


C IV (all) 


Mbh - 


L 


13795 


0.0257 


2.33E-03 


C IV (SDSS) 


Mbh - 


z 


11861 


0.0215 


5.97E-02 


C IV (SDSS) 


Mbh - 


L 


11861 


0.0255 


9.94E-03 


C IV (2QZ) 


Mbh - 


z 


1593 


0.0341 


5.71E-01 


C IV (2QZ) 


Mbh - 


L 


1593 


0.0303 


6.78E-01 


C IV (2SLAQ) 


Mbh - 


z 


307 


0.0642 


7.48E-01 


C IV (2SLAQ) 


Mbh - 


L 


307 


0.0572 


8.24E-01 



^ Number of quasars in each sample. 

^ Mg II velocities and BH masses corrected as described in Eqs. [2] 
and [3] 



Information from quasar broad emission lines 



5 



o - Hfi: 0.1<z<0.4 45.0<L„ ,<45.5 
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A<M>= 0.358 ( 2.14(j) 
N„= 12 
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- H(S: 0.4<z<0.7 45.0<Lb^|<45.5 
<M(j> = 8.301 +-0.009 rms=0,41 
<Mp>=8.334+-0.009 rms=0.39 
A<M>= -0.033 ( 2.62ct) 
N„=2109 




_ HI3: 0.4<z<0.7 45.5<Lb^|<46.0 
<Mq> = 8.562 + -0.006 rms=0.36 
<M^>=a.554+-0.007 rms=0.39 

" A<M>= 0.008 ( 0.93tj) 
N„=3623 




iog(M„H) (Mq) 



H^: 0.4<z<0.7 46.0<Lb_j|<46.5 
<Mo>=8.923+-0.015 rms=0.33 
<MB> = 8.a23H — 0.018 rms = 0.39 
A<M>= 0.101 ( 4.32(j) 
N„= 483 



log(M„J (MJ 



Figure 4. The distribution of black hole mass estimates in different redshift and bolometric luminosity intervals for H /3. We compare 
the original mass estimates (Mq; black solid lines) to the randomized estimates (A/r; red dashed lines). Mean masses are given for both, 
along with the difference in the mass and its significance. 



uncertainties on the virial estimators are typically 0.3- 
0.4 dex. Vestergaard & Peterson (2006) find a scatter 
of 0.43 and 0.33-0.43 dex for H /3 and C IV respectively, 
while McLure & Dunlop (2004) find a scatter of 0.33 dex 
for the comparison between Mg II and H /3 based mass 
estimates. This suggests that taking a mean velocity 
width to estimate a BH mass based simply on luminosity 
could be a useful approach, all be it, with some increase 
in the related error on BH mass. 

The distribution of Mbh vs. redshift is shown in Fig. 
[3] for the H /3, Mg II and C IV lines. The original dis- 
tributions are shown in black and the distribution after 
randomizing the emission line velocity widths are shown 
in red. In all three cases the original and randomized 
black hole masses appear essentially indistinguishable. 
To examine whether there is a quantitative difference we 
perform a 2-D Kolmogorov-Smirnov (KS) test (Peacock 
1983). We carry out these tests in both the Mbh — z 
and Mbh — L planes. To provide a robust estimate we 
generate 100 random realizations of the randomized BH 
masses and then take the median Dks and probability, 
although taking a single realization does not significantly 
change our conclusions. Table [1] contains the results of 
the 2-D KS analysis. For the H /3 line there is no sig- 
nificant difference between the original and randomized 
distributions of Mbh — z or A/bh — L (the null hypothesis 
that they are drawn from the same distribution is only 



rejected at the ~ 20% and ~ 10% level for Mbh — z or 
Mbh — L respectively). Applying this test to the Mg II 
mass estimates we find similar values of Dks i but due to 
the larger number of quasars the difference between the 
original and randomized BH masses is now highly signif- 
icant. Careful inspection of Fig. [3] does show that the 
randomized Mg II BH masses have a slightly broader dis- 
tribution at high redshift than the original BH masses. 
This could be due to intrinsic correlations of Mg II line 
width with redshift or luminosity (which indeed are seen 
by Fine et al. 2008). Alternatively, it could be driven 
by correlations of measurement errors with redshift or 
luminosity. Larger measurement errors will broaden the 
BH mass distribution, and if signal-to-noise is correlated 
with luminosity or redshift this will cause a difference be- 
tween our original and randomized BH masses. If we ap- 
ply corrections to the Mg II velocities as described in Sec. 
[21 then the significance of the difference in the distribu- 
tion of original and randomized BH masses is reduced to 
only 1-2(7. However, to be conservative we will consider 
the Mg II mass estimates without this correction for the 
remainder of the paper. For the C IV line the derived 
Dks values are again similar, and the significance of the 
difference is at the 2 (Mbh — z) to ia (Mbh ~ L) level. 

It might be expected that a combination of the Ed- 
dington limit and the bright end of the BH mass function 
might naturally cause a narrowness in the distribution of 
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_ Mgll: 0.4<z<0.7 45.0<Lg^|<45.5 
<M(,> = B.507+-0.007 rms=0.36 
<Mp>=8.582 + -q,006 rms=0.31 
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Mgll: 0.4<z<0.7 46.0<Lg^|<46.5 
<Mq>=9.1 19 + -0.017 rms=0.25 
<M„> = 9.124+-0.021 rms = 0.31 
A<M>= -0.005 ( 0.18cj) 



logCMg^) (Mq) 



Mgll: 0.7<z<1 .1 45.0<Lgjj|<45.5 
<M(j>=8.603+-0.007 rms=0.35 
<M„>=a.614+-q,00B rms=0.31 



A<M>= -0.011 i ).17ct) 
- N„=24ie 




Mg 



0.7<z<1.1 45.5<L„„,<46.0 



<Mq> = 8.870 + -0.004 rms=0.32 
■ <M^>=a.883 + -0.004 rms=0.31 
A<M>= -0.013 ( 2.20a) 
N„=5962 
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iog(M„H) (Mq) 



Mgll: 0.7<z<1.1 46.0<Lb^|<46.5 
<Mo>=9.1 35+-0.004 rms=0.26 
<Mp> = 9.135H — 0.005 rms = 0.31 
A<M>= 0.000 ( O.OOcj) 
N„=3430 




log(M„J (MJ 



Mgll: 1.1<z<1.5 45.0<Lb^|<45.5 
<M(,> = B.633+-0.01 5 rms=0.39 
<Mp> = 8.631 +-0.01 2 rms=0.30 
A<M>= 0.002 ( 0.12a) 
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A<M>= 0.008 ( 1.32a) 
N„=5056 




Mgll: 1.1<z<1.5 46.0<Lg^|<46.5 
<MQ>=g.1 9a + -0.003 rms=0.26 
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Figure 5. The distribution of black hole mass estimates in different redshift and bolometric luminosity intervals for Mg 11. We compare 
the original mass estimates (A/q; black solid lines) to the randomized estimates (Afj^; red dashed lines). Mean masses are given for both, 
along with the difference in the mass and its significance. 



velocity widths. To examine this we carry out the same 
analysis as above, but this time keeping each sample sep- 
arate. The results are listed in Table [TJ In the case of 
the 2SLAQ and 2QZ samples, there is no significant dif- 
ference between the original and randomized BH masses. 

Our second test is to measure the mean and RMS 
BH masses using each estimator in various redshift and 
bolometric luminosity bins. Again we calculate the val- 
ues for the randomized sample from 100 realizations, al- 
though we calculate the error on the mean assuming a 



single random sample. The results of this analysis are 
presented in Figs. IH [5] and [6] In all cases there is 
a good agreement between the original (Mq) and ran- 
domized (Mr) mass estimates. For H (3 the difference 
in mean mass is always less than 0.1 dex (except for 
the 0.1 < z < 0.4, 46 < LboI < 46.5 interval which 
only contains 12 quasars) and the median difference is 
~ 0.02 — 0.03 dex. Three z — Lbo\ intervals have differ- 
ences in the mean masses which are greater than 2cr. The 
comparison of Mg II BH masses in Fig. [5] shows similar 
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CIV: 1,5<z<1.9 45.5<Lg^|<46.0 
" <M(,> = 9.1 10+-0.012 rms=0.25 
<Mp>=9.128 + -0.01 0rrms=0.20 
A<M>= -0.018 ( I.Vdcj) 
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CIV: 1.5<z<K9 46.0<Lg^|<46.5 
<Mj,>=9.31 6 + -0.003 rms=0.21 
■ <M„>=9.324+-0.003 rms=0.20 
A<M>= -0.008 ( 1.83ct) 
N.=4023 , 
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CIV: 1.9<z<2.6 45.5<Lg^|<46.0 
<Mq>=9.009 + -0.020 rms=0.22 
<M„> = 9.106+-0.018 rms = 0.20 
A<M>= -0.097 ( 3.56cj) 
Nn= 121 

R, 



CIV: 1.9<z<2.6 46.0<Lb^|<46.5 
<M(j>=9.300+-0.004 rms=0.22 
<Mp>=9.331 +-0.004 rms=0.20 
A<M>= -0.031 ( 5.78ct) 
N„=3045 




CIV: 2.6<z<3.3 46.0<Lbjj|<46.5 
<Mq>=9.349 + -0.030 rms=0.28 
<Mp>=9.396 + -0.021 rms=0.20 
A<M>= -0.047 ( 1 .28a) 
N„= 88 



CIV: 3.3<z<4.0 46.0<Lgjj|<46.5 
<Mo>=9.460+-0.034 rms=0.24 
<Mp> = 9.395H — 0.028 rms = 0.20 
A<M>= 0.064 ( 1 .46(j) 
N„= 50 



log(M„H) (Mq) 



Figure 6. The distribution of black hole mass estimates in different redshift and bolometric luminosity intervals for C IV. We compare 
the original mass estimates (JV/q; black solid lines) to the randomized estimates (A/r; red dashed lines). Mean masses are given for both, 
along with the difference in the mass and its significance. 



good agreement. The largest difference is 0.075 dex, and 
the median absolute difference is only 0.01 dex. Again, 
there are several intervals where the difference between 
the mean masses is significant. The C IV masses (Fig. [5]) 
are consistent with the other emission lines. The maxi- 
mum difference is 0.1 dex and the median is 0.03 — 0.05 
dex. Two of the intervals have significantly (> 2a) dif- 
ferent means. Overall, while there are some significant 
differences in the mean Mbh, there is impressive agree- 
ment between the original and randomized BH masses. 
Such agreement has a number of implications, which we 
will discuss below. 

4. DISCUSSION 

The virial method makes use of the radius-luminosity 
relation derived from reverberation mapping, so that the 
only observables required are velocity width and lumi- 
nosity. Thus we have 



BH 



A X L" X ctblR' 



(4) 



where A is a normalizing constant, Lx is the monochro- 
matic luminosity at some wavelength A near the emission 
hue, a is ~ 0.5 (e.g. Bentz et al. 2009) and ctblr is the 
width of the emission line. There are multiple stages re- 
quired to make the connection between BH mass and the 
directly observable values. Starting with BH mass these 
are: 



1. Black hole mass (Mbh) and accretion rate 
(L/^Edd) give us bolometric luminosity, LboI- 

2. LboI combined with bolometric corrections and 
photometric errors give us L\, the first observable. 

3. LboI (or possibly iuv) defines the radius of the 
broad line region i?BLR- 

4. Mbh and i?BLRj together with the assumption of 
Keplerian orbits, gives us the velocity of the broad 
line gas, Vblr- 

5. Vblr together with non-virial motion, orientation 
effects and measurement errors give us ctblr,, the 
second observable. 

An important point to note is that, apart from uncer- 
tainty in the bolometric correction and photometric er- 
rors, scatter or errors in each of these steps will increase 
the scatter in ctblr at a given luminosity. The measured 
uncertainties in virial relations are typically 0.3 — 0.4 dex 
(Vestergaard & Peterson 2006), and these are only sta- 
tistical errors and do not include contributions from sys- 
tematic uncertainties in the reverberation mapping data, 
to which the virial relations are calibrated. The differ- 
ence between the mean BH masses measured in our orig- 
inal and randomized samples was at most 0.1 dex, and 



8 



Croom 



typically 0.02 dcx. This is much smaller than the uncer- 
tainties in the estimators themselves. Indeed to highlight 
this fact, we compare the SDSS Mg II BH masses calcu- 
lated by Shen et al. (2008) and Fine et al. (2008). Al- 
though both use the same virial relation calibration, the 
diflFerence in method (FWHM vs. IPV) results in a typ- 
ical difference of the mean BH mass (in narrow redshift 
and luminosity bins) of ^ 0.25 dex. This is an order of 
magnitude greater than the difference between our orig- 
inal and random estimates, although it is worth noting 
that the line measurement methods of these two works 
are somewhat different. 

The above results suggest that assuming a mean ve- 
locity width and an RMS scatter about that mean would 
enable us to estimate the BH mass distribution to the 
same precision as we currently have using the individ- 
ual measured line widths. This has a number of im- 
plications for the use of virial BH masses. One of the 
most common derived quantities is the accretion effi- 
ciency L/LEdd- As the velocity width adds little informa- 
tion, L/L-Edd L^~", where a ~ 0.5. In a flux limited 
sample the most distant objects will be the most lumi- 
nous and so will be found to have the highest efficiency. 

We know that quasar broad emission line properties do 
correlate with various observables. The best known case 
being the Baldwin (1977) effect which correlates equiv- 
alent width with luminosity. Fine et al. (2008) finds 
that although the mean width of the Mg II line does 
not depend on luminosity, the dispersion of line widths 
docs. There is a significantly broader distribution of line 
widths at faint luminosities. No doubt this contributes 
to the small but significant difference we do see between 
our original and randomized BH masses. In contrast. 
Fine et al. (2010) find no correlation between the dis- 
persion in C IV width and luminosity (once corrected 
for measurement errors). The low overall dispersion in 
line width found (< O.ldex at high luminosity) by both 
these works naturally infers that in a flux limited sam- 
ple quasar broad line widths cannot have a large im- 
pact on BH mass estimates. The small scatter and other 
observed properties (e.g. the presence of strong out- 
flows indicated by broad absorption lines) may infer a 
problem with our assumption that the broad line gas is 
virialized. However, there are a number of observations 
which suggest that virialization is at least partially cor- 
rect. The most compelling being the correlation between 
lag and line width (e.g. Peterson et al. 2004). Compar- 
isons between broad line reverberation mapping masses 
and bulge velocity dispersion shows a correlation anal- 
ogous to that observed in quiescent galaxies (Onken et 
al. 2004). Added to this, there is agreement between 
dynamical mass and reverberation masses for the few 
galaxies where measurements have been possible (Davies 
et al. 2006; Onken et al. 2007; Hicks & Malkan 2007). 

The small impact of velocity width on BH masses is due 
to the limited dynamic range of velocity widths. Under 
assumption that the virial method is giving us real infor- 
mation concerning BH mass, the small range in velocity 
widths must be in part due to joint constraints from a 
steep BH mass function at high masses and the Edding- 
ton limit. Fine et al. (2008) examine this and find that 
some bias is possible, but that selection effects cannot 
on their own produce the narrowness of the velocity dis- 



tribution for the most luminous quasars. A more robust 
examination would be to extend the approach of Schulze 
& Wisotzki (2010) to the larger SDSS, 2QZ and 2SLAQ 
samples, although this is outside the scope of the present 
paper. However, we do find that there is no significant 
difference between original and randomized BH masses 
when we test the fainter 2QZ and 2SLAQ samples sepa- 
rately. The dynamic range in luminosity sampled by the 
combined SDSS, 2QZ and 2SLAQ sample is over 6 mags 
for the Mg II line (i.e. a factor of ~ 250 in luminosity). 
Given that L/L^dd ~ L^'^ this corresponds to a factor of 
^ 16 change in L/L^dd- It is rather remarkable that over 
this range there is no change the mean velocity width in 
the Mg II line. 

Clearly, if a sub-sample is selected on the basis of line 
width, then randomizing those velocities with the rest of 
the population will have significant impact on the derived 
BH masses. For example, narrow-line Seyfert Is (NLSls) 
are rare in any flux limited sample and havc^ bcnni consid- 
ered to have abnormally low BH masses. Such random- 
ization as we carry out will clearly change dramatically 
the mass estimates of these objects. However, recent 
work has suggested that NLSls do not have low masses 
but that their narrow lines are caused by orientation (De- 
car li et al. 2008) and/or radiation pressure (Marconi et 
al. 2008). Such effects should increase the scatter in ve- 
locity width for a given black hole mass. However, the 
impact of this is limited by the small observed scatter as 
a function of luminosity. 

Specific quasar sub-samples have been shown to have 
different virial BH mass distributions. One particular 
example is radio-loud or radio-detected quasars. Jarvis 
& McLure (2006) show that a sample of radio quiet 
quasars has narrower broad lines than a matched radio- 
loud sample. This demonstrates that the difference be- 
tween radio loud and radio quiet quasars is not due to 
selection biases. However, such biases do exist. The 
fraction of radio-detected quasars increases with lumi- 
nosity (e.g. Jiang et al. 2007), as one would expect in a 
flux limited sample, so any radio-loud sample will have 
higher estimates for BH mass even before accounting for 
the intrinsic property differences. 

A positive outcome of the weak impact of quasar 
broad-line width on black hole mass estimates is that 
a reasonable mass estimate can be derived without even 
measuring the width of quasar emission lines! With im- 
proved multi-band photometry (e.g. including near- and 
mid-IR data, Richards et al. 2009) high quality, faint, 
photometric samples of quasars can be obtained, with 
reasonable photometric redshifts. It would then be pos- 
sible to derive a quasar BH mass function without ob- 
taining spectra. However, inferences from such analyses 
would be limited. Basing BH masses simply on observed 
luminosity does not allow us to investigate accretion rate 
variations. 

A clear route to improve our current position it to 
attempt better calibration of the virial mass relations, 
and in particular determine how other parameters influ- 
ence line width (e.g. orientation, radiation pressure). A 
key question, which has yet to be answered, is whether 
there is any reason for the coincidence that accretion ef- 
ficiency scales with luminosity (as L/L^dd ~ L^'^) such 
that quasar broad line widths are constant with luminos- 
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